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Summary. We introduce an enhanced model based on the generalized centrifugal 
force model. Furthermore, the desired direction of pedestrians is investigated. A 
new approach leaning on the well-known concept of static and dynamic floor-fields 
| in cellular automata is presented. Numerical results of the model are presented and 

O f compared with empirical data. 

Force-based models try to describe the dynamics of pedestrians as reaction 
to forces acting on each single pedestrian. Basically two kinds of forces can 
be distinguished: 

• driving forces designed to drive pedestrians to a desired direction with a 
desired speed. 

• repulsive forces which are responsible for preserving the volume exclusion 
of pedestrians. 

Since the introduction of force-based models [1] many works were dedicated 
to investigations of the repulsive forces and finding new and better forms [2-6] . 
These efforts for improving the form of the repulsive force is understandable, 
since the interactions between pedestrians dominate the dynamics, especially 
at high densities. Surprisingly, not much work has been done on the influence 
of the specific form of the driving force which is expected to dominate the 
behavior at low or intermediate densities. 
The standard form of the driving force is 

H drv =m*i^, (1) 
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with a relaxation time r and a desired velocity . Although this expression 
is simple, it is not clear how to choose the desired direction 

IK II 

in a given situation and only very few works were concerned with modeling the 
desired direction (2). In [7] an Ansatz with directing lines was introduced to 
steer pedestrians around 90° and 180° corners. Gloor ct al. [8] used a path- 
oriented approach to model the desired direction of agents on given hiking 
paths. 

In [9] Moussai'd et al. have formulated the determination of the desired 
direction in form of a minimization problem. 

It should be mentioned that the directing problem we discuss here, i.e. the 
determination of the desired direction for each pedestrian, is conceptually dif- 
ferent from the classical routing problem. In [10] an algorithm for generating 
automatically a navigation graph in complex buildings in combination with 
directing lines at corners was proposed. Another algorithm for way finding in 
buildings was proposed in [11]. Recently a further development of the notion 
of the "quickest path" using a non- iterative method to estimate the desired di- 
rection in the social force model (SFM) was introduced [12]. The main concern 
in this class of problems is how to define and connect intermediate targets, in 
order to facilitate the evacuation of pedestrians. By contrast, in the directing 
problem the existence of such intermediate targets is in general assumed. 

In this work we introduce enhancements of the generalized centrifugal 
force model (GCFM) and investigate on their basis the modeling of the desired 
direction (2). For the sake of demonstration we test our model in two different 
geometries: a bottleneck and a corner. 



1 The model 

In this section we give a brief overview of the GCFM and its definition. Fur- 
thermore we introduce an effective modification of the pedestrian-wall inter- 
actions that simplifies the definition of the repulsive force. 

1.1 Pedestrian-pedestrian repulsive interactions 

Introducing the vector connecting the positions of pedestrians i and j, 

Rij Rj Rit &ij ) ; (3) 



Rij 



the repulsive force in the GCFM reads 
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F{j ^ — Tfiikij cij , (4) 

with rrii = 1 the mass of i and the effective distance between pedestrian i and 

3, _^ 

d ij =11 II -n(wi) - rj(vj) , (5) 

and the polar radius rj of pedestrian i. 

The relative velocity is defined such that slower pedestrians are less 
affected by the presence of faster pedestrians in front of them: 

1, 



2 -[(vi - Vj) -e tj + \{ Vi -Vj)-e tJ \ 



(vi - vj) ■ eij if (vi - vj) ■ e i} > 
otherwise. 



(6) 



The parameter 



fCj, ■ 



2 Vi 

T7> 



(vi ■ e7j)/ || Vi || if W • > & || W ||^ 
otherwise, 



(7) 



reduces the effective range of the repulsive force to the angle of vision. Through 
the coefficient kij the strength of the repulsive force depends on the angle: it 
is maximal when pedestrian j is in the direction of motion of pedestrian i and 
minimal when the angle between j and i is bigger than 90° . 



1.2 Wall-pedestrian repulsive interactions 

In the GCFM the interactions between pedestrians and walls are modeled by a 
force similar to the pedestrian-pedestrian repulsive force. A wall is represented 
by three point masses acting on pedestrians within a certain range. From 
a computational point of view this analogy exhibits an overhead since the 
repulsive force between a pedestrian and a wall is calculated three times. 

We now make use of the "distance of closest approach" as defined in [6] 
to formulate the repulsive force between a pedestrian i and a wall w as 

f£™p = rf || ? || k iw b iw , (8) 

with 

^="(i-^0-(i-^i). 

where I is the distance of closest approach between an ellipse and a line, r 
is the polar radius determined by the nearest point on the line to the center 




Fig. 1: Illustration of distances used in the definition of the wall-pedestrian 
repulsive force (8). 



of the ellipse i (Fig. 1). HQ is the Heaviside step function, ki w is defined in 

Eq.(7), || Vf* || is the desired speed of i and rf is a parameter to control the 
strength of the force. 

The repulsive force (8) is a contact force that is different from zero if the 
effective distance of the center of the ellipse to the segment line is non-positive. 
For the simulations in this paper we set the strength of the repulsive forces 
as rj — 0.2 and rf = 5. 



2 Influence of the desired direction 

In this section we study the effects of the desired direction on the dynamics 
of a system by measuring the outflow from a bottleneck with different widths. 
See Fig. 2 for the simulation set-up. Four different methods for setting the 
direction of the desired velocity are introduced and discussed. Finally, simu- 
lation results will be compared. 




Fig. 2: Scenario set-up: Pedestrians move from a holding area (shaded region) 
through the bottleneck (I = 2.8m, h = 4.5m, b = 4m and w variable). 
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2.1 Strategy 1: Directing towards the middle of the exit 

The first strategy is probably the most obvious one. Herein, the desired di- 
rection ef for pedestrian i is permanently directed towards a reference point 
that exactly lies on the middle of the exit. In some situations it happens that 
pedestrians can not get to the chosen reference point without colliding with 
walls. To avoid this and to make sure that all pedestrians can "see" the middle 
of the exit the reference point ei is shifted by half the minimal shoulder length 
b m m = 0.2 m (Fig. 3). Pedestrians that pass to the right of the reference point 
ei head towards e2- 



Fig. 3: Strategy 1: All pedestrians are directed exactly towards the reference 
points ei and e^. 

2.2 Strategy 2: Enhanced directing towards the middle of the exit 

This is a modification of strategy 1. Pedestrian are still directed to the 
shifted reference point e\. However, from a certain position pedestrians can see 
through the bottleneck the second reference point e-x- In this case e\ is ignored 
and the desired direction is set to be parallel to the line e~\e\- Since, pedestri- 
ans that are inside the bottleneck can always see ei the desired direction is 
kept parallel to e\e\. 

Here again the reference points and the delimiting range of the bottleneck 
is shifted in x- and y-direction by 6 max (Fig. 4). 
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Fig. 4: Strategy 2: Depending on their position pedestrians adapt their di- 
rection. In the range where the exit of the bottleneck is visible (marked by 
dashed lines) the direction is longitudinal. Outside this area they are directed 
towards the middle of the bottleneck. 

2.3 Strategy 3: Directing towards the nearest point on the exit 

Another possibility to choose the desired direction ef is to define a line I in 
front of the exit and take at each time the nearest point from the pedestrian 
i to I (Fig. 5). In comparison with strategy 2, pedestrians that are not in the 
range where the point e 2 is not visible choose one of the end points of the line 
/. In strategy 2 this would be the middle of I. 




Fig. 5: Strategy 3: Directing towards the nearest point on the exit. Molnar 
published in [13] a very similar strategy. The only difference is the placement 
of the line, which is away from the corner by 6 ma x- 

2.4 Strategy 4: Guiding line segments 

Without loss of generality we introduce the main idea of strategy 4 with 
help of the previous bottleneck. Unlike the previous strategies this strategy is 
applicable to all geometries with corners even if the exit point is not visible. 
We recall that in strategy 3 a line in front of the bottleneck was defined. The 
nearest point from each pedestrian to this line was set to define the desired 
direction. As a generalization we make use in strategy 4 of three different lines 
to "smoothen" merging in front of the bottleneck (Fig. 6). 
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Fig. 6: Strategy 4: Guiding line segments in front of the bottleneck. For each 
corner a set of three line segments is generated. The length of all directing 
lines is equal to 3.5 m. 

The blue line set (down the dashed line segment) is considered by pedes- 
trians in the lower half and the red line set by pedestrians in the upper half 
of the bottleneck. For a pedestrian i at position pi we define the angle 



with j such that 0j = min{#i, 02, #3}. As in strategy 3 the direction lines are 
shifted in x- and y-direction by 6 m i n . 

2.5 Numerical results 

In the previous section we have proposed different methods (called strategies) 

for choosing the desired direction ef . To compare these strategies we have 
performed simulations for a bottleneck using the same set of parameters for 
the GCFM. For each strategy only the width of the bottleneck was varied 
from 1 m to 2.4 m. 

On the basis of a quantitative analysis the importance of the choice of 
strategy for the observed behavior can be estimated. In the following, for 
each strategy we measure the flow through bottlenecks of varying width w. 
The flow is measured directly at the entrance of the bottleneck according to 




(10) 



with hj the nearest point of the line j to the pedestrian i. 
The next direction is then chosen as 



~J) _ Pihj 




(11) 



J = 



N At -l 
At 



(12) 



with iV/i t = 60 pedestrians and At the time necessary that all pedestrians 
pass the measurement line. 
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In Fig. 7 the resulting flow for all four strategies is presented. 




w [m] 

Fig. 7: Flow through a bottleneck with different widths. Simulation results 
with different strategics for the desired direction of pedestrians in comparison 
with empirical data from [14]. The experiments were conducted with 180 
persons. 

The flow for strategy 1 saturates independently of the width. This was 
expected since pedestrians do not use the whole width of the bottleneck and 
keep indeed oriented to the middle. The picture changes for strategies 2-4, 
where the effective width of the bottleneck is clearly larger. Strategy 2 shows a 
better usage of the middle widths (< 1.8m). Here, the slight blocking near the 
corners, that emerges from strategy 3 is particularly disadvantageous. Strategy 
4 produces higher flows for widths up to 2 m. Up 2 m the flow stagnates. The 
main observations are: 

• The choice of the strategy for the desired direction influences considerably 
the outcome of the simulation. 

• An inconsiderate choice of strategy, in that case strategy 1, can lead to 
large variations from experimental results. 

• In contrast to strategy 1, strategies 2-4 show better usage of the bottle- 
neck width and lead to higher flow values. 
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3 An application: Motion around a corner 

Basically, force-based models are functional only in areas, where the exit is 
constantly visible by all pedestrians. Obviously, this can not always be guar- 
anteed which is a problem since a proper initialization of the desired direction 

ef for each pedestrian is not possible. In order to overcome this problem one 
has to introduce "virtual" exits. This was showcased previously with strategy 
4. 

In this section, we introduce enhancements of strategy 4 and study their 
impact on the movement time, i.e. the time until all pedestrians have left 
the simulation set up. For simplicity we consider the movement of TV = 100 
pedestrians in a 90°-corner-like corridor. 

The basis of our enhancements is the following observation: Given a guid- 
ing line I, the desired direction of a pedestrian i is determined in dependence 
of its position and the nearest point to I. This choice neglects two important 
factors: 

1. The perception of space: Individuals try to minimize, when possible, their 
path to the exit. In our example, pedestrians would take a point near the 
corner as goal and not the nearest point on the guiding line. Depending on 
the starting position of pedestrians, this can be far away from the corner 
and much longer than the shortest path to the exit. 

2. The dynamical and collective influence of pedestrians: In the presence of 
other pedestrians and depending on the magnitude of the local density, 
the nature of the "quickest path" [15] changes dynamically and differs in 
most cases from the "shortest path" to the exit. 

We therefore adopt a concept similar to ideas introduced in [16] which are 
well established and widely used in cellular automata models [17, 18]. 

At a time step t a pedestrian i heads towards a point on the line which 
minimizes the distance to the inner point of the corner l l {t). This is a natural 
territorial effect which leads to the shortest path to the exit. If all pedestrians 
try to take the shortest path, large jams will be observed right at the inner 
point of the corner. If, however, the collective influence of pedestrians domi- 
nates the choice of the desired direction, pedestrians will choose their desired 
direction to be orthogonal to the guiding lines and thus make better use of 
the whole directing line. 

For this reason we include a dynamical factor that depends mainly on 
previous decisions taken by other pedestrians: 



p % (t) =cxp(-fc d -occ* i(t)) , 



(13) 



where 



occJ el (t) = 



(14) 



in ax 



is a measure of the occupation of the line. n l is the cardinality of the set 
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A l = {P \ j eB t kP <P} 



and n\ 



max 



is the cardinality of the set 



Bi = {j e [l, N] I i± j k ^i-E^> 0}. 



(15) 



Bi is the set of all relevant neighbors of i, that influences its desired direction 
by means of a contribution to occj. ol (t) (14). For the scenario depicted in Fig. 8 
the set Bi for i (red ellipse) contains only one pedestrian j (bold ellipse). 



Fig. 8: How to get around the corner? Pedestrian i that is heading toward 
the first guiding line, considers the positions of its neighboring pedestrians as 
well as its initial position to decide whether or not to head closer to the edge 
of the corner. 

Large values of occ* ol (£) imply small values of p l (t). As a consequence 
pedestrians prefer not to change their desired direction closer to the edge of 
the corner. 

Finally, the update rule of the distance l*(t) is given by: 



p l {t) € [0, 1] gives the rate of change from the initial "guess" of pedestrian i. 
For p*(t) = the desired direction of i stays orthogonal to the guiding line, 
while p l (t) = 1 displays the case where i's desired direction is directed to 
the edge of the corner A. In the next section we study the influence of the 
parameter kd on the dynamics of pedestrians. For the second and third line 
we set kd — and vary it only for the first line. 




a 

J A 




(16) 
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4 Analysis of the sensitivity parameter 

To understand the impact of the collective influence of pedestrians on the 
chosen target point for each pedestrian i, we study the time evolution of the 
relative length for different values of kd- The relative length is defined as 

iL(t) = (17) 

Wax 

where l max is the length of the guiding line. 

Fig. 9a shows the probability distribution of the relative length for kd = 0. 
Pedestrians are mainly heading towards A and the full length of the directing 
line is rarely used. 



0.0 0.2 



0.6 0.8 



MR 

0.0 0.2 



0.6 0.8 



(a) k d = 



(b) k d = 10 



Fig. 9: Movement of N = 100 pedestrians around a corner with different values 
of kd- The width of the corridor is w — 3 m. 



The situation changes considerably for kd = 10. Fig. 9b shows that the 
distribution of the length is more balanced which indicates that pedestrians 
make better use of the directing lines. 

To showcase the impact of collective influence of pedestrians on the desired 
direction, we show in Fig. 10 the variation of the movement time in dependence 
of k d . 
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Fig. 10: Movement time for a simulation with N — 100 pedestrians around a 
corner for different values of kd- 

A qualitative comparison shown in Fig. 11 confirms the above-mentioned 
quantitative analysis. 




Fig. 11: Screen shot of a simulation with 100 pedestrians, kd — (top) and 
kd = 10 (bottom). 

For kd = a jam forms immediately before the corner as indicated by the 
large number of slowly moving pedestrians (red ellipses). This results from 
a strong competition between the pedestrian to pass close to the edge A of 
the corner. In contrast, for kd = 10 pedestrians move quicker since they make 
optimal use of the guiding line. 
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5 Summary 

Wc have developed a strategy to determine the desired direction ef for each 
pedestrian i. This method is rather general and can be used in each geom- 
etry characterized by the existence of corners, e.g. bottlenecks (2 corners), 
T- Junction (2 corners). In analogy to CA models we introduced and tested a 
factor to model the static and dynamic interactions of pedestrians with the 
geometry. 

Our work was based on an enhanced version of the GCFM [6]. The en- 
hancements use a considerable simplification of the repulsive forces acting 
on pedestrians from walls. Furthermore, we addressed an important issue in 
force-based models, namely the choice of the desired direction of pedestrians. 
Several strategies were implemented and compared with empirical data. This 
comparative investigation showed that the outcome of a simulation depends 
strongly on the chosen direction of the desired direction of pedestrians. Fi- 
nally, we introduced a new mechanism to direct pedestrians in 90°-corners by 
means of directing lines. The main concept of this strategy base on the well- 
known concept of dynamical floor-field. For further works, the parameter kd 
that expresses the tendency of pedestrians to take the shortest path (or not) 
should be varied individually as the geometrical and dynamical conception of 
pedestrians differ. 
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